Field theory of bi- and tetracritical points: Statics 
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We calculate the static critical behavior of systems of 0(nu)(B0(nj_) symmetry by renormalization 
group method within the minimal subtraction scheme in two loop order. Summation methods lead 
to fixed points describing multicritical behavior. Their stability boarder lines in the space of order 
parameter components ny and n± and spatial dimension d are calculated. The essential features 
obtained already in two loop order for the interesting case of an antiferromagnet in a magnetic field 
(nil = 1, n± = 2) are the stability of the biconical fixed point and the neighborhood of the stability 
border lines to the other fixed points leading to very small transient exponents. We are also able to 
calculate the flow of static couplings, which allows to consider the attraction region. Depending on 
the nonuniversal background parameters the existence of different multicritical behavior (bicritical 
or tetracritical) is possible including a triple point. 



PACS numbers: 05.50.+q, 64.60.ae 



I. INTRODUCTION 

Antiferromagnets in an external magnetic field show a 
variety of phase diagrams depending on the interaction 
terms present in the spin Hamiltonian^ . The spin in- 
teraction may be isotropic, anisotropic with an easy axis 
and/or single-ion anisotropy terms, where the anisotropy 
is in the direction of the external magnetic field. The 
phase diagram of such models exhibit a multicritical 
point, where several transition lines meet. 

At a bicritical point three phases - an antiferromagnetic 
phase, a spin flop phase and the paramagnetic phase - are 
in coexistence. The phase transition lines to the param- 
agnetic phase are second order transition lines, whereas 
the transition line between the spin flop and the anti- 
ferromagnetic phase is of first order. At the tetracritical 
point four phases - an antiferromagnetic phase, a spin flop 
phase, an intermediate or mixed phase and the paramag- 
netic phase - are in coexistence. All transition lines are 
of second order in this case. 

A field theoretic description of these models starts 
with a static functional for an rt-component field $ of 
0(n») © 0(n±) symmetry (nil + n± = n) leading to 
different multicritical behavior connected with the sta- 
ble fixed point (FP) found in the renormalization group 
treatment^^. Bicritical behavior has been connected 
with the stability of the well known isotropic Heisen- 
berg fixed point of 0(n\\ + symmetry, whereas tetra- 
criticality has been connected with a fixed point of 
0{n\\) © 0(n±) symmetry, which might be either the so 
called biconical FP or the decoupling FP. In the last FP 
the parallel and the perpendicular components of the or- 
der parameter (OP) are asymptotically decoupled. 

The important questions which theory should give an 
answer to is, (i) which of these FPs is the stable one 
in a three- or two-dimensional system, and (ii) what are 



the differences in the critical behavior at the multicritical 
point? These questions have been risen and considered 
in one loop order—, where the Heisenberg FP turns out 



to be the stable one in d = 3 for the case 
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1 and 



n±_ = 2, but this picture is changed in higher loop order. 
In a five loop order e = 4 — d expansion it has been found 
that the biconical FP is the stable one^. It also has been 
found that the differences between the exponents at the 
different multicritical points are much smaller than in the 
one loop order calculation. 

Physical examples where such multicritical behavior 
has been found are the anisotropic antiferromagnets^ 
(with the magnetic field in the hard direction) like 7,8 
GdAL0 3 ano£ MnF 2 , as well as±£ MnCl 2 4D 2 or-ii 
Mn 2 AS4 (A=Si or Ge). Other examples with a single 
ion anisotropy might be layered cuprate antiferromag- 
nets like (Ca,La)i4Cu 2 404i. Besides the examples with 
rt|| = 1 and n± = 2 one might consider other cases: 
1 and n± 
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= 1 when additional anisotropies are 
present as i n 12 i 13 NiCl 2 or high-T c superconductors rep- 
resenting a system with — 2 (corresponding to the 
superconductor OP) and n± = 3 (corresponding to the 
antiferromagnetic OP). 

Quite recently the possible types of phase diagrams in 
the magnetic field - temperature plane of d = 3 of uniax- 
ially anisotropic antiferromagnets have been studied by 
Monte Carlo simulation o 14 i 15 . For n\\ = 1 and n± = 2 
a phase diagram with a bicritical point has been found 
in agreement with earlier simulations^, but contrary to 
the results of renormalization group theory in higher loop 
orders^. 

A general picture is obtained when one considers a 
generalized model with an n component order parame- 
ter (OP), which splits into n|| parallel OP components 
and nj_ perpendicular OP components and quartic inter- 
action terms of 0(n\\) © 0(n±) symmetry. Both paral- 
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lei and perpendicular OP components become critical at 
the multicritical point. In the n||-n^-space regions of dif- 
ferent types of multicriticality exist touching each other 
at stability border lines ('phase border lines') where the 
fixed points change their stability (such a picture of the 
different stability regions might be called a 'phase dia- 
gram'). In addition to the stability of a fixed point we 
want to mention that one has to consider also the at- 
traction regions of a fixed point to answer the question 
wether one can reach the stable fixed point. In order 
to discuss the attraction regions one has to consider the 
flow of the couplings from the nonuniversal initial (back- 
ground) values. 

We therefore reconsider the critical behavior of systems 
with O(nii) © 0{n±_) symmetry. Being interested in crit- 
icality of three dimensional systems, we will work within 
the minimal subtraction scheme and evaluate the results 
at fixed dimension d = For the universal proper- 
ties (as asymptotic critical exponents and marginal di- 
mensions) it turns out that already the two loop cal- 
culations refined by resummation are in good quantita- 
tive agreement with previous resummed higher order e- 
expansion results. 5 However, contrary to previous calcu- 
lations, the technique we use gives a possibility to an- 
alyze non-universal effective critical behavior which is 
manifested in a broader temperature interval near the 
(multi)critical point. Such calculations are out-of-reach 
the £-expansion and will be performed below on the base 
of analysis of the renormalization group flow. 

The paper is organized as follows: starting from the 
static functional (Sec. [TTJ) we introduce the renormaliza- 
tion in Sec. IIIII and calculate the field theoretic functions 
in Sec. IIVI The perturbative expansions being asymp- 
totic, we apply in Sec. |V]the resummation technique to 
restore their convergence and to extract numerical values 
of the fixed points of the renormalization group trans- 
formations and their stability. We discuss the stability 
border lines between the fixed points and show that they 
are shifted considerably compared to the one loop cal- 
culation. As a result for the isotropic antifcrromagnet 
in a magnetic field represented by the point (1,2) in the 
n||-n^-space the biconical fixed point is stable predicting 
tetracritical behavior if the fixed point is reached from 
the background. The very neighborhood of the stability 
border lines is characterized by very small transient ex- 
ponents. Moreover, looking at the attraction regions a 
surface in the space of the fourth order couplings exists 
above which no finite fixed point can be reached. This 
indicates the possibility of a scenario mentioned already 
earlie r 5 ' 18 ? 19 where the multicritical point is a triple point 
and the second order lines separating the paramagnetic 
phase from the ordered phases contain a tricritical point. 
In section fVTl the critical exponents are defined and their 
asymptotic values are calculated for the physically inter- 
esting case ri|| = 1, n±_ — 2. The flow equations and 
effective exponents are discussed in Sec. IVIII leading to 
our final conclusions and outlook in Sec. IVIIIl In the 
appendices we discuss the perturbative expansion for the 



vertex functions (Appendix A) and explain the resumma- 
tion procedure exploited in our calculations (Appendix 
B). 



II. STATIC FUNCTIONAL 

The critical behavior of an isotropic system (0(n) sym- 
metry) with short range interaction is determined by the 
static functional 

f ( 1 ^ ^ l " 

I i=l 

+ J(Vo-0o) 2 j, (1) 

which is known as the Ginzburg-Landau-Wilson(GLW)- 
functional. The order parameter (f>o = (j)o(x) is assumed 
to be a n-component real vector. The symbol • denotes 
the scalar product between vectors, r is proportional to 
the temperature distance to the critical point and u is 
the fourth order coupling in which perturbation expan- 
sion is usually performed. Systems represented by such 
a static functional have been extensively studied in the 
last decades with different renormalization procedures, 
and the corresponding critical exponents and amplitude 
ratios are known up to high loop orders (see e.g.—). 

In order to describe bicritical behavior the n- 
dimensional space of the order parameter components 
will be divided into two subspaces with dimensions nj_ 
and n|| with the property rtj_ + ?i|| = n in the following. 
Correspondingly the order parameter separates into 



f_L0 

P\\o 



(2) 



where </>j_o is the n_L-dimensional order parameter of the 
n^-subspace, and tf>»Q is the rz.ii -dimensional order param- 
eter of the ny-subspace. Performing the separation in the 
GLW-functional (p} one obtains 



n 



Bi : 



d xl -r±4>j_ 



1 _ _ 

i=i 

nil 

1 -» -> 1 ' _ _ 

2^11^110 ' Ho + 2 ^2 ^^\\o ' v ^i|o 

8=1 

(0.LO • 0J-oJ + 4}- (^>||0 • Ho) 



This functional contains three fourth order couplings 
{u} — {u±,u x ,u\\} and instead of one parameter f as 
in ([T]), in ([3]) two different parameters r± and rn appear 
referring to different temperature distances. 
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The decomposition in parallel and perpendicular OP 
components allows to describe the critical behavior at the 
meeting point of two critical lines: (i) the line where r± 
becomes zero and the n^-dimcnsional components </>j_o 
are the OP, and (ii) the line where f|| becomes zero and 

the n\\ -dimensional components </>|| are the OP. At the 
meeting point both quadratic terms become zero and 
both components of <po have to be taken into account. 
The critical behavior of this multicritical point, has been 
described already in one loop order— and three differ- 
ent types of multicritical behavior have been found; (i) 
one described by the well known isotropic n component 
Heisenberg fixed point where all fourth order couplings 
are equal, (ii) one described by a decoupling fixed point, 
which consists of a combination of two n± and nil com- 
ponent isotropic Heisenberg fixed points of two decou- 
pled systems and (iii) a new type of fixed point called 
the biconical fixed point. Which of these fixed point is 
the stable one depends on the number of components 
n± and n« and the dimension d of space. The scaling 
properties depend on the symmetry of stable fixed point 
When the 0(n) symmetry in the OP space is broken to 
0(ni|) © 0(n±) symmetry, then also the spatial correla- 
tions are different for the two OP subspaces. 

III. RENORMALIZATION 



would be possible to introduce a 3 x 3-matrix for the 
renormalization of the three fourth order couplings {u} 
as performed ir£. The resulting /3-functions are the same 
anyway. Thus it remains a matter of taste if one uses a 
scalar or a matrix renormalization for the three fourth 
order couplings. The situation changes if one considers 
(f) 2 - insertions. The vertex function r^ 2 ' 1 ) splits up into 
four functions tf'^ ± , tf'lL fftj| and ffn.'V (for the 
notation see Appendix [A"|). A consistent renormalization 
of the 4> 2 - insertions is only possible by introducing 

with the renormalization matrix 

»<■ = ( £ £ ) • « 10 > 

The zeroth order of the perturbation expansion only ap- 
pears in the diagonal elements Zu = 1 + 0({u}), while 
the off-diagonal elements Yij = 0({u}) start with the one 
loop order. With (j4|) and ([9]) an arbitrary vertex function 
(|A9|) renormalizes as 

r (N,L) _ 7 l/2 yl/2 

X H ( Z <p 2 )pm ■ ■ ■ { Z ^)pLiL^ { a 1 - L) aN -t 1 --- lL ( n ) 



The procedure used to obtain the vertex functions ap- 
propriate for the renormalization is described in more 
details in the Appendix [Al From the general structure of 
the two point vertex functions presented therein follows 
that the order parameter functions </>j_o and 0|| O in the 
subspaces may be renormahzed by the scalar renormal- 
ization factors 



= Z 



V2, 



= z 



V2; 



(4) 



The above relations and the definitions (|A7|) and (|A8|) im- 
ply that the correlation lengths, £|| and do not renor- 
malize. They constitute together with the wave vector 
modulus k the independent lengths of the system. The 
fourth order couplings may also be renormahzed by scalar 
renormalization factors 



uj_ = k 6 ' Z^Z u ^u^A d 1 , 
u x = K e Z^\Z^Z u ^u x A^ x , 
H = ^ZT^Z^A- 1 , 



with the geometrical factor 



2/ V 2/ (27r) d 



(5) 
(6) 
(7) 



(8) 



analogously to DohmSi. In ([5]) T(x) denotes the Eu- 
ler Gamma-function and fid is the surface of the d- 
dimensional unit sphere. We want to remark that it 



The indices aj, (3j and ij are running over _L and ||. 
The renormalization introduced in the above relations 
removes the poles from all vertex functions except the 
functions T.^^ 2 which have to be considered separately. 
Thus equation (fTTj) is valid for all vertex functions except 



.(0,2) 

;/3i/3 2 



The multiplicative renormalization in (|lip does 
not remove all singularities at d — 4 in these function. 
In order to remove the remaining poles an additional ad- 
ditive renormalization is necessary. The multiplicative 
renormalization leads to functions 



.(0,2) 



f(0,2) 

2*2 ;i!i2 



(12) 



A finite function can be obtained by subtracting the sin- 
gular part [T^:f i02 } s from (H]). Thus we may introduce 

r§& = «' V (rgg* - ir&Js) (13) 

as renormahzed functions. The additive renormalizations 
are then defined by 



Ap lP2 {{u}) = ^A- d \T^ ih ] s 



(14) 



The three functions A\\\\ and A^\\ = A\\±_ in 



([14)1 represent the extension of the function A{u) in the 
isotropic case 2 - and may be written as the symmetric 
matrix 



A({u}) = 



A± 
A 



A x 
A,, 



(15) 
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Within statics it may also be convenient to work with 
the temperature dependent vertex functions without in- 
troducing the correlation length, as described in the Ap- 
pendix at stage (|A6[) . This allows to avoid functions with 
</> 2 -insertions except for the specific heat. In this case a 
renormalization for the temperature distances Arj_ and 
Afii has to be introduced as performed in2i. In order 
to obtain renormalization factors which do not contain 
ratios of the temperature distances a matrix renormal- 
ization 



Afj_ 

'II 



^ I )=Z^-Zr-Ar (16) 



with the matrices 



Z r = 



Z r± Y r 



Zrb — 



z+ x 

Zj,, 



(17) 



has to be used. The renormalization factors in the matrix 
above are obtained by collecting the £-poles proportional 
to Arj_ and Arii in the two vertex functions f and 
f ^' 0) . The matrices in (JT7J) are related to the renormal- 
ization matrix Z ^ for the 2 -insertions by 



— (Z i X Z r 



zlzi} 



(18) 



which represents the matrix counterpart of the scalar re- 
lation in the isotropic case. The superscript T denotes 
the transposed matrix. 



IV. C - AND f3 - FUNCTIONS 

From the scalar renormalization factors Z$ a . (a, =_L 
|) the (- functions 



&«,({«}) = 



d In k 



(19) 



are derived. The K-derivatives, also in the following defi- 
nitions, always are taken at fixed unrenormalized param- 
eters. In two loop order we obtain the (^-functions 



of the matrix are 

nj_ + 2 



[C/ 



u± ( 1 - —tii. ) - ^ul (23) 
6 V 12 / 72 x v ' 



ia i "II fi M x 



(24) 
(25) 



= ^ll(l-^ll)-^ (26) 

The non diagonal elements (f2"4")) and ([23)) are proportional 
to h x and differ only in a prefactor n±_ and nil respec- 
tively. Thus we may write 

[C^]i2=nx« x C(M) = ^[C^] 2 i (27) 



with a function 



(28) 



in two loop order. Relation (fT8")) between Z^ and Z r 
implies 



C<*2 — Cr — Co 



(29) 



where the diagonal matrix C,^ = diag(^_ L , Cd>n ) has been 
introduced. From the additive renormalization (|15p the 
function 



[B^i]ij({u}) = K £ ^\Z \p2\jn\Z 



tf> 2 \jp 



n , p 



x^- K - £ ^ ! 1 ]„ J [Z f! 1 ] pm [A] Im ({ H }) (30) 



l.m 



can be introduced which is the extension of the scalar 
function (u) in the isotropic case (see Ref.— for defi- 
nitions). Calculating (f3"0")) in two loop order we obtain 

B^(W)=(f m)+0({« 2 }) (31) 
The /3-functions of the four couplings are defined as 



72 



72 x 



72 



(20) 
(21) 



From the matrix Z \2 of the renormalization of the 4> 2 - 
insertions it is convenient to introduce the ^-matrix 



Mii(W) = -E(^[^ 2 ]«)[^ 



(22) 



A matrix C, r following from Z r , which has been defined 
in (fTT)) . also can be introduced analogous to (|22p by re- 
placing Z 02 with Z r . In two loop order the components 



Pu a ({u}) = K—^j— 

an 



(32) 



with a =_L, ||, x. In two loop order the explicit expres- 
sions of the /3-functions are 

(n±+8) 2 nu „ 

DO 

- (3 ^ + 14) ,j - ^ Ul _ul - X* ; (33) 
12 - 1 36 x 9 x ' v ' 

(ni+2) (n,| +2) 
Pu x = su x H liiiix H 5 u x it|| 



6 



2 (ni + n|| + 16) 3 (ni + 2) 2 



72 



5 



+2) 



uiun 



5(n ± + 2) 2 
72 U ^ 



5(n,|+2) 2 

-U x Ui| , 



72 



/3«„ = -euii 



("II + ■ 
6 



1 2 , n± 2 
+ — « 



X 



(3tt.ii + 14) 5n± 2 ni o 
-U\\ unut M y 

|| 3 g II X g X 



12 



(34) 



The flow equations of the fourth order couplings u a are 



(35) where I is the flow parameter. 
I 



(36) 



FP 




«*x 


u| 
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TABLE I: Fixed points and stability exponents of the 0(n«) © 0(n±) model. 
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V. FIXED POINTS AND THEIR STABILITY 
('PHASE DIAGRAM') 

The FPs of the flow equations Eq. (|3"6"|) are given by the 
solutions of system of equations: 



(37) 



At the one loop order, Eq. (|3"T)) defines eight fixed points. 
Six of them are real and two are complex apart from the 
region where nn and n± are small (see table H]). They all 
have the property to be proportional to e = 4 — d. To 
proceed with higher loop approximations one can make 
use of different calculation schemes, either performing 
an e-expansion or solving the flow equations directly for 
fixed d = 3 (i.e. for e = A particular feature of the 
e-expansion is that an increase of the order of approxima- 
tion does not lead to an increase of the number of FPs^ 
Once the FPs are found in the first order of e, the next 
orders of the expansion only give the next order contri- 
butions to the first order values of the FPs but do not 
lead to the appearance of new FPs. On contrary, when 
one directly solves non-linear flow equations for fixed d 
in higher-loop order, more and more fixed points may 
appear in addition since the order of the polynomials to 
be solved for the fixed points increases. Moreover, the 
physical FPs found in e-expansion may disappear when 
one naively solves the flow equations for fixed space di- 
mension d. It is well established by now that the ex- 
pansions involved for the field theoretic renormalization 
group functions are asymptotic at bes1i2£ and one has to 
use appropriate resummation procedures to get reliable 



results on their basis. Let us note however that whereas 
in the e-expansion the resummation procedure is needed 
only to precise the FP values, whereas in the fixed d ap- 
proach it is the resummation that allows to judge about 
presence of the FP at all and comparative analysis of the 
two approaches allows to judge about the FP picture on 
a sound basis^. Below we will make use of these two 
approaches. We will work within the two- loop approx- 
imation and show that the two loop e-expansion is not 
sufficient even when one uses resummation in contrast to 
the fixed dimension procedure. 

Defining the FP picture we are looking for answers to 
the questions: (i) is the critical behavior described by 
a certain FP with corresponding asymptotic and uni- 
versal critical exponents? (ii) is the system very near 
a stability border line and slow transients leading to an 
effective and nonuniversal critical behavior? In such a 
case then the whole information contained in the nonlin- 
ear flow equations is necessary and available only in the 
second method presented here. Furthermore, in answer- 
ing these questions we demonstrate that the physically 
relevant features are obtained already using the two loop 
approximation. 



A. Results in two loop e-expansion 



Starting from the Eqs. (f33|) - (f35|) one gets six real FPs, 
as shown in the table [U Four fixed points correspond 
to the decoupled effective Hamiltonians of 0(n±) and 
0(n\\) models. Besides the Gaussian FP G(u* ± — u* = 
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uf, = 0) these are the H(n±) and H(n\\) FPs with (u* ± — 
u H ^\ u* x = u| = 0) and (uj { = u n{ - n w\ u* x = u* ± = 0), 
correspondingly, as well as the decoupling FP 2?(w^ = 
u n( - n± \ u* = 0, ut = u H ^). Here and below by u n ^ 
we denote the Heisenberg FP of the 0(n)-symmetrical 
model. Two more FPs correspond to the non-zero value 
of the coupling u x . Following the nomenclature of Ref, 4 
we call them the isotropic Heisenberg and biconical FPs, 
H(n± + and B correspondingly (see table UJ). In the 

minimal subtraction RG scheme, the value of 

u H(n) ig 

currently known with the e 5 order^ From any of the 
/3-functions (f33|) . (|35|) one recovers familiar e 2 result: 



,n(n) _ 



6e 18(3n+14)e 2 



(n + 8) 3 



(38) 



Expressions for the coordinates of the FP B are too cum- 
bersome to be given in a compact form for general n« and 
n± . Below we list the non zero values for the FP accord- 
ing to Tab. U for the physically important case n\\ = 1, 
n± = 2 we are interested in. Note, that for this case the 
two FPs U\ and U2 attain complex values: 

u* {1) = 0.66667e + 0.41975e 2 , 
u™ (2) = 0.60000e + 0.36000e 2 , 
u* {3} = 0.54545e + 0.31104e 2 , 
tt| = 0.40496e + 0.38691e 2 , 
ul = 0.50569e + 0.31610e 2 , 
u® = 0.69059e + 0.29926e 2 , 

u" 1 = (0.41521 - i 0.32936)e + (0.21259 - i 0.21159)e 2 , 
(0.20213 + i 0.12400)e + (0.15182 + i 0.7087)e 2 , 
(0.98646 + i 0.12302)e + (0.62025 + i 0.06112)e 2 



The bad convergence of the e-expansion is already shown 
in the two loop order values of the biconical FP B in 
d = 3. This FP does not fulfill the criterion 



A* = ulu* A 



,*2 



> 



(39) 



for describing tetracritical behavior as in zero loop order 
(see Eq. (5.12) in Ref4 and the discussion in Sec. IVIII 
below) . 

From the structure of the /^-functions one can derive 
exact values of some of the stability exponents, as shown 
in the table |U There, by we denote usual stability 

exponent of the 0(n) model. The rest of the exponents 
are defined in the appropriate FPs by: 



H(n A 



= dp Ux /du x \ H(n±) , 
= d/3 Ux /du x \ H{ni[) , 



= dp Ux /du x \v- 



(40) 
(41) 
(42) 



To find the stability exponents in the FPs TC(n± + n\\) 
and B one should solve the appropriate secular equation. 

One can see from the table that three FPs Q, Tt[n±), 
and W(nii) are never stable for d < 4. However as we 



will see below, the stability of the other three FPs for 
d < 4 depends on the values of n± and . As far as two 
stability exponents of the decoupling FP V are always 
positive, ^("J./'"!' > 0, the stability of V is defined 
by the sign of the exponent uJ§ . Therefore, the equation 
for the marginal dimensions line in the nj_ — ri||-plane 
reads: 



d/3 Ux /du > 



IP) 



0. 



(43) 



Equation (|4*3")> defines a curve n^(n\\) (or, equivalently 
nV(nx)) that borders a region of n±,ri\\ values where 
the FP T> is stable. Substituting the second-order result 
into the function (fM)) we get from (|4*3")l : 



ra? (ny) 



2(16 -?i||) 



48e 



Eq. (|4"4")l can be inverted and one gets: 



nV(n±) 



2(16 -n x ) 



48e 



n ± + 2 n ± + 2 



(44) 



(45) 



The first term in (|44|) . (I45|) coincides with the first order 
result of Ref. 4 whereas the second term gives the second- 
order contribution, which again demonstrates the weak- 
ness of the second order e-expansion, since the shift to 
smaller values of Ti||(nx) in d = 3 is drastically overesti- 
mated leading to instability for small values of n\\ and/or 

Tlx- 

Note, that the stability properties of the FP D can be 
evaluated on the base of exact scaling arguments.— At 
this FP the coupling term u x 0x^ji nas a scaling of the 
product of two energy-like operators, the latter having 
scaling dimensions (l — a n± )/v n± and (1 — a n ,,)/v nt , cor- 
respondingly (with a„, v n being the heat capacity and 
the correlation length critical exponents in 0(n) univer- 
sality class). In turn, this leads to the following formula 
for the RG dimension y Uy . of the combined operator^ 2 ^ 



(46) 



With available five-loop e-expansion results for the expo- 
nents of 0(n) theory 2 ^ the marginal dimensions rt?(rty) 
(nV(n±)) can be estimated. 

The first order result for n?(ri||) is shown in the figure 
[T]by a dashed line. FP T> is stable for the values of n±, n\\ 
above this line. Crossing the line, one gets into the region 
where the biconical FP B acquires stability. With further 
change of nx,^|| its coordinates do change as well and 
for certain value of n?(n||) this FP coincides with the 
Heisenberg FP H(n^ + n±). Then it looses its stability 
and further the FP Ti(n^ + n±) is stable. Therefore, the 
marginal dimension line n?(n||) (or n^(n±) equivalently) 
can be defined from the any of the conditions: 

(3 UJ _ (ux,uj_,u_l) = Pu x (u x ,u x ,u x ) = 

/3 U|| (M||,U||,U||) = 0. (47) 
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FP 


Ml 


w* 


u\ 


Wi 


U>2 


LJ3 


g 











-1 


-1 


-1 


W(2) 


1.141 








0.581 


-0.461 


-1 


W(l) 








1.315 


-1 


-0.552 


0.565 


V 


1.141 





1.315 


0.581 


-0.014 


0.566 


W(3) 


1.002 


1.002 


1.002 


0.597 


0.407 


-0.036 




1.128 


0.301 


1.287 


0.583 


0.554 


0.01 



TABLE II: Fixed points and stability exponents of the O(l) © 0(2) model obtained by the Pade-Borel resummation within 
two loops. Biconical FP B is stable. 



As far as in the Heisenberg FP TL the RG functions de- 
pend on the sum of field dimensionalities n±_ + n\\, the 
resulting marginal dimension curve will be of the form 
n±(n>\\) = const — n\\. Substituting the FP TL coordi- 
nates into expressions for the /3-functions (f3"3")) - (f3"5")) we 
get the second order expression for the marginal dimen- 
sion. By it we recover the two first terms of the third 
order result quoted in Ref4 (and obtained from Ref. 27 ): 

n±{n\{) = -n|| +4-2e + c x e 2 , (48) 

with c x = ^(6C(3) - 1) = 2.5885. Or, equivalently: 

nf(n±) = -n± +4- 2e + c x e 2 . (49) 

In Fig. [T]the results obtained in the e-expansion in first 
loop order are shown. Whereas the 7iZ?-stability border- 
lines lead to an acceptable result, but bad convergence, 
the B2?-stability borderline show unphysical features in 
the e-expansion. In second order in e the stability bor- 
derline for positive values of n± lies at negative values 
of 7i|| meaning that the decoupling fixed point is stable 
in the whole region shown. The second order e expan- 
sion results do not lead to reliable results being estimated 
naively. Therefore, below we will reanalyze the RG func- 
tions by resuming them in two loops directly for s = 1. 

The FPs where the parallel and perpendicular system 
decouples (u* = 0) need some comments. The renor- 
malization group procedure used here assumes that the 
multicritical system is described by one diverging length 
scale and therefore by one correlation length £ and one 
corresponding critical exponent v. This does not hold for 
decoupled systems where two length scales and there- 
fore two correlation lengths, £m and £j_ with two different 
asymptotic exponents and v±_ are present. Thus the 
usual scaling laws with one length scale break down (see 
the remarks in Ref. 4 and Ref. 28 ). 

B. Results from resummation of the two loop field 
theoretic functions 

Let us pass now to another way of analysis of the RG 
functions (|55)) - (|55l) . Namely, in the spirit of the fixed 
dimension RG approach 1 ^ let us consider the FP equa- 
tions ([57)) directly at fixed d — 3. The RG series being 




1 2 3 4 5 



n par 



FIG. 1: Regions of different static bicritical behavior in the 
nil — n^-plane (e = 4 — d = 1), which are defined by the 
stable FP (Heisenberg FP TL, biconical FP B and decoupling 
FP T>) . Shown are the 7iZ?-stability borderline (dashed lines) 
and 6D-stability borderlines (solid lines), in one loop order 
(thin lines) and resummed two loop order (thick lines). The 
dots indicate low integer values for OP components. One sees 
that the borderlines are drastically shifted to smaller values 
of OP components. Thus in the case n\ = 1 and n± — 2 FP 
B (connected with tetracriticality) is stable in two loop order 
contrary to the one loop order calculation where the FP TL 
(connected with bicriticality) is stable. 



divergent, we present the /3-functions (|3"3")) - (|3"5"|) in a form 
of resolvent series and resum them by the Pade-Borel re- 
summation technique as explained in the Appendix [B] 
Let us note that such a representation preserves symme- 
try properties of the functions. Now, the FP coordinates 
as well as the marginal dimension lines are evaluated nu- 
merically. The lines n^(n||), n^(ny) are shown in Fig. [T] 
by solid lines. The FP coordinates for n± = l,rt|| = 2 
and the stability exponents are given in Table [TXJ 

From the resummation at n\\ =1 and different n±_ we 
can follow the changes in the FP values of the fourth or- 
der couplings at the biconical FP (see Fig. They 
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start with equal values corresponding to the isotropic 
Heisenberg FP at the borderline value n — 1.6. The cross 
coupling between the parallel and perpendicular compo- 
nents decreases to zero at the stability borderline to the 
decoupling FP at n±(n« = 1) = 2.17, whereas the other 
couplings slightly increase. 



1,4 






i 




1,2 










1,0 










0,8 










0,6 










0,4 


■ - ■ -u 

par 








0,2 


U 

perp 

U* 

X 








0,0 








1,7 1,8 


1,9 


2,0 


2,1 



'perp 



FIG. 2: Dependence of fixed point values of the static cou- 
plings u_t_, uii , u x (in Borel summed two loop order) on n± for 
ri|| = 1 in the region of stability of the bicritical fixed point 
in d = 3. 



VI. IDENTIFICATION OF THE CRITICAL 
EXPONENTS 

The connection between critical exponents and the £- 
functions can be obtained from the solution of the renor- 
malization group equations for the vertex functions. 



A. Anomalous dimensions 7/ and r/± 

Considering the vertex functions (|A9p at the bicritical 
point £J 2 = = the renormalization group equation 
for the fc-dependent functions read 

d „ d 



N 



du a ' 2 . 

a=_L,||,X 3 = 1 

x r£:?? aN (k,{u},K) = o. (50) 



Solving the equation with the method of characteristic 
equations leads to 



6n 



dx 
x 



x r 



3 = 1 J1 

™(^M0})(5i) 



where 5^ = N + d — \Nd is the naive dimension of the 
vertex functions. The couplings {u(l)} are determined by 
the flow equations (|36j) . For N — 2 the vertex functions 
represents the ^-dependent inverse susceptibilities in the 
two subspaces. Ea. (|5Tj) reduces to 



r£'51(fc,W,«) = («0 a «q) 



dx 



i x 
k 



?™(^ v Ml)}) (52) 



with the indices (not to confuse with the specific heat ex- 
ponents below) cx,i equal to _L or ||. In the asymptotic re- 
gion the couplings {u(l)} have nearly reached their fixed 
point values {u*}. The (^-function C<fi a ({ u ( x )}) m the 
exponential factor reduces to the constant ({m*}) = 
C • Expression ([52")) reduces to 



r&2(*.M.*) = * a i a+<J -*^(^K} 



(53) 



In order to obtain a finite amplitude function f L^a- we 
choose the matching condition 



k 
kI 



= 1 



(54) 



which determines the flow parameter l(k). Inserting into 
(1531 leads to 



1 K 



2 +C 



fg'S (LK}) 



(55) 



Therefore the asymptotic behavior of the inverse suscep- 
tibilities Xam = ^ 



k^***! ~ k 2 -^ . (56) 
Thus we may identify the two anomalous dimensions 

1°, = ~Q ai (57) 

with oii =_L, || . 

B. Exponents of the susceptibilities 7 and 7j_ 

In order to obtain the exponents 7|| and 7^ and their 
relation to the exponent of the correlation length we 
have to consider the temperature dependent vertex func- 
tions and the corresponding renormalization group equa- 
tions. There are two methods to include the temper- 
ature dependence into the renormalization group equa- 
tions. The first one is to constitute renormalization group 
equations for the vertex functions (|A6|) and derive from 
the renormalization of the temperature distance Ar a 
relation between the exponents 7y and 7j_ and the corre- 
sponding function £ r . The second one is to constitute a 
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renormalization group equation for the vertex functions 
(|A9p including the temperature dependence by an ex- 
pansion in 2 -insertions. This would lead to a relation 
of the exponents 71 and 7j_ to the function £j 2 . As a 
consequence of relation (fT8|) one obtains the scaling laws 
between the exponents 7M and 7_l and m and rj± with in- 
cluding only one exponent v for both correlation length 
(see below). In the following we will consider the first 
method. The renormalization group equations for the 
vertex functions (IA6I) at k = read 



d 



£ a. 

a=_L,||,x 



d 



d 



dA 



3=1 



du a 

r { a N i : Q) aN ({Ar},{u},K) = (58) 



where relation (I29|) already has been used. The matrix 
C^2 can be diagonalized by the transformation 

diag(C + ,C-)=^ 1 Cj 2 P (59) 

where and C- are the eigenvalues of C,d> 2 , while the ma- 
trix P is determined by the corresponding eigenvectors. 
The matrix has the form 

1 / V± n^u x C \ 



t T - 

SiP- — 



6 \ n±u x C V\\ 



(60) 



The functions V± = V±({u}), V\\ = Vj|({u}) and C* = 
C({u}) may contain all orders of loop expansion. The 
two loop expressions have been given in [)23 p -(|28 p . The 
eigenvalues of this matrix read 



C± = i(£±w) 

where we have introduced the square root 
W = 

A = Vj_ — Vj| . 



v 'A 2 + 4n_Ln||U 2 C 2 
and the parameters 

£ = V± + Vj| , / 



(61) 



(62) 



The corresponding eigenvalues constitute the transfor- 
mation matrix 



1 

2n±Ux C 
A+W 



2ny u x C 
A+W 
1 



(64) 



where the first and the second column are the eigenvec- 
tors to C+ and C— • The explicit appearance of the pa- 
rameters (|62|) and (|63|) in the matrix (|64|) may differ by 



using the relation 

(A + W)(A-W) = -4n ± n||M 2 C 2 . (65) 
For special cases the eigenvalues and the eigenvectors 
may simplify considerably. 

i) In the case of the decoupled system, where u x = 0, 
a decay into two isolated isotropic subsystems with or- 
der parameter dimensions n+ and occur. The ma- 
trix (f64|) reduces to the unit matrix and the eigenvalues 

C+ = V±/6 = and C- = Vf|/6 = (p } are the corre- 

sponding isotropic ^-functions of the two subsystems. 

ii) In the case of the isotropic Heisenberg system, i.e. 
u 1 = My =ii» — u, the matrix (1641) reduces to 



(66) 



1 



which is independent on u also in higher order perturba- 
tion expansion. The first eigenvector cf^_ = (1, 1), which 
corresponds to £ + , points in a 45 degree direction from 
the bicritical point in the Arj_ -Arn -plane, while the sec- 
ond eigenvector <f^_ — (— nn/nj_, 1) lies in some sense tan- 
gentially to it. The eigenvalues in two loop order are 

n+ + njj +2 / 5 \ 
C+ = g « I 1 - ^ u J = , (67) 



u / n_i_ + nil + 10 
C = 3 { 24 u 



(68) 



from which one can see that £ + reproduces in this case 
the corresponding ^-function of an isotropic n± + nii- 
component Heisenberg system. The Eigenvalues and the 
direction of the Eigenvectors imply that obviously de- 
fines an exponent is, while the other Eigenvalue £_ is 
connected with a crossover exponent. 

Diagonalizing and introducing transformed tem- 
perature distances 



r± 



= P' 1 A : 



(69) 



(63) in Eq. (|58p we obtain the renormalization group equation 



d \ _ d d d 

K x~ + — ^^+ r +a — — 

ok * — ' ou a or + OT- 

a=X,||,x 



1 W 



3=1 



rL 1 / .'°L JV (r + ,r_,{ U }, K ) = 0.(70) 



The solution of this equation is 
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N A 



T<£:°i N (r + ,r-,{u}, k) = («M N exp [± £ f 

.7 = 1 Jl 



dx 

x 



-Lai- 



^(N,Q) f r+(l) MO 
Qn V(k0 2 ' (kI) 2 



,{«(!)} 



(71) 



r 



where the transformed temperature distances fulfill the 
flow equations 



^ = r ±C± (W). 



(72) 



Considering the solution ([7T|l in the asymptotic region 
{it} = {u*} we obtain for N = 2 



T^(r + ,r-,{u},K) = K 2 l 



2,2+fL 



x 



(2,0) ( MO MO r n \ f73) 



As discussed above, r + is the temperature distance to 
the bicritical point, thus the matching condition 



MO 



= i 



(74) 



defines the flow parameter l(r + ) as a function of temper- 
ature. In the asymptotic region the solution of Eq. (j?2"|) 

is 



MO = M c ± 

From (f?4"]) we obtain therefore 



/ 



(75) 



(76) 



Inserting into (173|) leads to the asymptotic expression 



rg^)(r +I r_ I W, K ) = /, 2 (^) 7 " < 

x«(i l7 ^,K}) 



(77) 



of the inverse order parameter susceptibility and we may 
identify the exponent r ) ai and the crossover exponent <f> 
as 



2 + G 



7a; 



<t>* 



2-a 



2-C- 
2-Ct 



(78) 



Introducing ?y Qi from (|5T[) into the first relation one can 
write 



7a, = ^ — = - »7aJ 



(79) 



2-0- 

which obviously leads to the exponent 

M = 2 - C* = M (80) 



of the correlation length. Although two exponents r\ ai 
(cti =_L, ||), and as a consequence two exponents j ai , 
have been obtained, only one exponent v describes the 
multicritical behavior. This reflects the fact that the 
anisotropy is present only in the order parameter com- 
ponent space, but not in the coordinate space. Only a 
single diverging length scale is present in the system 
in our case, for a discussion when two length scales are 
present see Reh 2 ^. For the discussion in section IVIDI it 
is convenient to define the exponent 



2-C* 



(81) 



quite analogous to ([BP) . The crossover exponent in (|75|) 
can then be written as 



(82) 



C. Exponent of the specific heat a 

Within the Ginzburg-Landau- Wilson model the spe- 
cific heat is proportional to the (/> 2 -0 2 -correlation func- 
tion, which is the negative vertex function r' ' 2 ). There- 
fore we have to consider the solutions of the renormaliza- 
tion group equations for rf^g a defined in p3[) in order 
to obtain a theoretical expression for the exponent a. 
Considering the functions rf2fg 2 as a symmetric matrix 



-1(0,2) _ 



r (0,2) p(0,2) 

1 ±_L A _L|| 

r (0,2) p(0,2) 

Mi ^ n n 



(83) 



the renormalization group equation reads 



a=±, || , X 



du a 



dAf 



,(0,2) 



,(0,2) 



+ r(°- 2 > • M = b 4 



(84) 



where the vertex function is taken at k = 0, i.e. r*- ' 2 -* = 
T^ ' 2 ' ({Ar}, {it}, k). The term — e is the naive dimension 
of the vertex function and it appears explicitly because 
r^ 2 -* has been introduced as dimensionless quantity in 
(fT3l) . Introducing the diagonalized (-functions (|59j) and 
the transformed temperature distances (fBT)]) , Eq. ([54]) can 
be rewritten as 



d v-^ n ° d j. d 

: ^~+ M P**n h M+^ M r -^ — 

ok i ou a or + or- 

a=_L,||,x 
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r^ 2) + 2dia g (c + ,C-)-r 



(0,2) _ 
± — 



B 



(±) 



The transformed vertex functions are 



.(0,2) 



r (0,2) p(0,2) 
p(0,2) r (0,2) 



Quite analogously is B^ 2 = P B^P. Applying the 

(0,2) _ 



T r (0,2), 



pip 



(85) 



(86) 



matrix ([66)1 from the isotropic case one obtains V 



isotropic Heisenberg model. Thus in the present case 
the specific heat is obviously proportional to r^°V^ . The 
corresponding renormalization group equation is 



d_ 

8 k 



a=_L, || , X 



d , d d 

7) vC >+ r +J) h< ^- r -^ — 

ou a or + or^ 



r 



(0,2) 



2r 



(0,2) 



rjj?, , which is the specific heat in the with the solution 



J 



-2C+ - e 



r (0,2) _ R (++) 
1 ++ - 



(87) 



r 



dx 



^- B { + + \{u{l)})e 



J'^gr (2C+-e) 



Using (|74|) - ([76]) the expression reduces in the asymptotic 
region to 



J 



T^\r + ,r-,{u},K) 



Vk 2 



^(0,2) 



1- 



7- 2 



( r+ / K 2)l/* 



(++) 



2c; 



(89) 



r 



where a is given by 

2-a 



a 



2-a 



2-dis. 



(90) 



In the last equality we have used Eq. ([50)l . Thus the hy- 
perscaling relation has been derived from the renormal- 
ization group equation for the specific heat. 



D. Numerical estimates for asymptotic critical 
exponents and scaling 

The results for marginal dimensions presented above 
in section [V] give an evidence that for the physically rel- 
evant case nil = 1, n± = 2 the biconical FP B is stable. 
Substituting the e-expansion for the biconical FP B we 
recover the results presented in Rcf. 5 . Although the sec- 
ond order of perturbation theory is known as an optimal 
truncation for the e-expansion, in our case the values for 
the exponents are not reliable. Resummation of the e- 
expansion in this order makes not much sense. 

Therefore we choose another way to estimate the nu- 
merical values of the exponents making use of the fixed 
d expansions. There, we proceed as follows. First, 



based on the two- loop expressions (|23|) - ([26)) for the ma- 
trix [C^a]*i({u}) P2"|) we find its eigenvalues as expansions 
in renormalized couplings. In turn, the exponents v+, i/_ 
are expressed in terms of the corresponding eigenvalues 
as series in renormalized couplings as well. Correspond- 
ing series are found also for the magnetic susceptibility 
exponents. The series are then resummed (as described 
in the Appendix [B| and evaluated at the FPs. Numeri- 
cal results of the exponents v+, i/_, 7| | , an d 7j_ obtained 
by such evaluation are given in table IIIII in the (stable) 
biconical FP B and, for comparison, in the (unstable) 
Heisenberg FP Ti. The expressions for the exponents rj± 
and m are too short to be resummed. Therefore their 
numerical values are found from the familiar scaling re- 
lations m = 2 — 711/1/4., r/± = 2 — ~f±/v + using resummed 
values of 7s and vs. For the decoupling FP one has to 
modify the scaling relations according to the statements 
above and has to use m = 2 — 711 /f_, rjx = 2 — jx/v+- 
Note that within the accuracy of calculations results for 
the exponents that correspond to parallel and perpendic- 
ular fields do not differ (the difference shows up within 
the fourth digit). The overall agrement with the 5 loop e- 
expansion is very good especially for the stable biconical 
FP. 
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Reference 


FP 


'/-L 


''II 


/-L 


Til 
'II 


7/ i 






(X 


this paper 


B 


0.037 


0.037 


1.366 


1.366 


0.696 


0.692 


1.144 9 


-0.088 




11(3) 


0.040 


0.040 


1.411 


1.411 


0.720 


0.564 


1.275 9 


-0.160 


Ref. 2 


B 








1.222 


1.222 


0.611 


0.503 


1.176 


0.167 


Ref^ 


W(3) 








1.227 


1.227 


0.611 


0.505 


1.136 


0.167 


Ref£ 


B 


0.037(5) 


0.037(5) 


1.37(7) 


1.37(7) 


0.70(3) 


0.56(3) 


1.25(1) 


-0.10(9) 


Ref^a 


W(3) 


0.0375(45) 


0.0375(45) 


1.382(9) 


1.382(9) 


0.7045(55) 


0.559(17) 


1.259(23) 


-0.114(17) 



9 Pole in the Pade approximant is present (see Appendix [EJ. 



TABLE III: Critical exponents of the O(l) © 0(2) model obtained by resummation of the two-loop RG series at fixed d = 3 
in different FPs (first two rows of the table). Our data is compared with the results of Ref.— (first order e-expansion), and 
Refs.— ^ (resummed fifth order e-expansion). Numbers, shown in italic were obtained via familiar scaling relations. 



VII. FLOW EQUATIONS AND EFFECTIVE 
EXPONENTS 

The resummation of the /3-function has the big advan- 
tage to find the fixed point values and asymptotic expo- 
nents but in addition the flow of the couplings from their 
background values to their FP values. The flow is impor- 
tant for the crossover behavior but also for the discussion 
of the asymptotic multicritical behavior of physical repre- 
sentatives of such systems with different the background 
values of the couplings. Their location in the attraction 
region of a FP defines the critical behavior.. 

A. Flow of fourth order couplings 

Fig. [3] shows the flow lines in the space of the coupling 
parameters for different initial conditions calculated from 
the resummed /^-functions. Mean field theory states a 
criterion- which has to be fulfilled by the fourth order 
couplings for the existence of a tetracritical point; it reads 

A = u\\u± -u\ > . (91) 

The flow equations show that A is not an invariant sur- 
face of the flow. However it contains several separatrices 
with the corresponding fixed points: the decoupling FPs 
with u* = ( Q, H(l), H(2), V) and the Heisenberg 
FP H(3). Therefore in the region shown in Fig. 1 the 
mean field condition A = represents quite well the sur- 
face (called MFS) separating the attraction region of the 
biconical FP B and the region of run away flows. 

Thus one can draw the following conclusions. A system 
with initial conditions A < lies outside the attraction 
region of the stable FB and its runaway flow indicates 
that a first order transition is expected. One would con- 
clude that the multicritical point is a triple point. Sys- 
tems with initial conditions where roughly A = would 
flow to the Heisenberg FP 7i(3) indicating that the mul- 
ticritical point is a bicritical point. Then finally if the 
initial conditions are such that A > the biconical FP 
B is reached and the multicritical point is a tetracriti- 
cal point. In this way the three scenarios sketched in 
Ref£ are realized. The important and open point is to 




0,0 



FIG. 3: Resummed flow for different initial conditions. The 
unstable FPs are shown as filled spheres, the stable biconical 
FP as filled cube. In order to show the crossover the ini- 
tial values of couplings are chosen accordingly. Due to the 
small transient exponent present the stable fixed point is not 
reached for the small value of the flow parameter chosen. The 
fixed points are connected by separatrices defining the surface 
which encloses the attraction region it is slightly different from 
the MFS-surface shown (for further details see text). The ini- 
tial values of flow 6 are outside the attraction regions. The 
flow parameter is changed in the interval —40 < InZ < 0. 

connect the initial conditions of the field theoretic flow 
equations to the interaction parameters in the appropri- 
ate spin Hamiltonian containing the anisotropic interac- 
tions defining the antiferromagnetic system. 



B. Effective exponents 

Having available the solutions of the flow equations one 
can define the effective exponents evaluating the field the- 
oretic ^-functions at the values of the couplings given by 
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the flow according to the definitions of the exponents in 
Eqs. ([57]) and (|50|) . We substitute the couplings obtained 
from the resummed flow equations into the resummed £- 
functions appearing in the expressions for exponents. In 
this way the effective exponents become functions of the 
flow parameter I, e.g. 



-1 



(I) = 2-C+(«y(0,«±(0,«x(0) = ^7/(0 
eff {l) = 2-C-(«||(0,ttx(0.«x(0) (92) 

ff follows from Eq. 



The effective crossover exponent 
dH as 



<t>eff{l) 



e//(0 



e//(0 



(93) 



Results for the exponents are shown in the Figs. EJ [5] 
and [5] 




FIG. 6: Effective exponents 
flows shown in Fig. [3] 



for the initial conditions of 



1,5,6 




In I 



FIG. 4: Effective exponent 711 (solid lines) and 7x (dashed 
lines) for the initial conditions of flows shown in Fig. [3] 



1,5 




In I 



FIG. 5: Effective exponents v for the initial conditions of 
flows shown in Fig. [3] 



The two exponents of the parallel and perpendicu- 
lar OP susceptibilities are almost equal even in the 
nonasymptotic region but might be quite different from 
the asymptotic value. Especially due to the slow tran- 
sient present the values for the flows 1-3 are smaller than 
the expected asymptotic values. The same holds for the 
effective exponent of the correlation length v e ff- The 
flow 6 which does not reach a finite fixed point lies within 
the expected values. The formal decrease of the effective 
exponents to unphysical negative values are due to the 
flow to infinite values of the couplings. 

VIII. CONCLUSION AND OUTLOOK 

We have shown that the two loop order perturba- 
tion theory together with appropriate resummation tech- 
niques is sufficient to calculate the multicritical behav- 
ior appearing in systems with O(ny) © 0(n±) symmetry 
at fixed dimension. The advantage of such a procedure 
lies in the accessibility of the corresponding flow equa- 
tions, which allow a discussion of attraction regions and 
effective (crossover) critical behavior. We confirmed the 
shift of the one loop stability borderlines, with the conse- 
quence that the multicritical behavior for the case n\\ = 1 
and n± = 2 is characterized by the stable biconical fixed 
point and not by to Heisenberg fixed point. The dis- 
cussion of the attraction region of this fixed point leads 
to the possibility of different phase diagrams depending 
on the nonuniversal initial parameters entering the flow 
equations. 

In a next step the results will be used to recon- 
sider dynamical critical behavior of the pure relaxational 
dynamics^ of these systems. For the bicritical dynam- 
ics of an antiferromagnet in an external magnetic field 
an extension of the above mentioned model is necessary 
since besides the OP conserved densities (CD) have to 
be taken into account. In statics these couplings only 
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appear up to second order terms and can be integrated 
out. In dynamics they lead to a coupling of the two 
dynamic equations. A complete description of the dy- 
namical multicritical behavior for has not been given in 
two loop order so far. Finally mode coupling terms have 
to be taken into account up to two loop order. 
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Project No. P19583-N20. 



APPENDIX A: VERTEX FUNCTIONS AND 
PERTURBATION EXPANSION 

With the static functional the vertex functions 



({*},{«},< 



(Al) 



can be calculated in a definite loop order by collecting 
all one-particle-irreducible contributions. N and L are 
the total numbers of <j>- and <f> 2 - insertions and the in- 
dices oti and f3j indicate if the corresponding insertion is 
of type _!_ or ||. {r} = r±, fii and {u} = u±, u x , un act as 
placeholder for the two temperature distances and three 
fourth order couplings. Between the lower critical dimen- 
sion d u = 2 and the upper critical dimension d c = 4 the 
vertex functions (|A1[) also contain singularities at d = 3. 
They have their origin in a non analytical shift of the crit- 
ical temperature as function of the four point couplings 
(for more details see the third refeerence in Refiii). n 
order to remove them, two parameters f± c and rii c are 
introduced. They are determined by 



f ( L 2 i 0) ({0},{u},0)=r Xc , 

-(2,0), 



r^({0},{«},0) 



(A2) 
(A3) 



This defines to functions r± c ({u}) and f^ c ({u}). Intro- 
ducing 



A rj_ = f± - r± c ({u}) , 
Af|| = r\\ - f\\ c ({u}) 



(A4) 
(A5) 



and rewriting the expressions for (|A1[) leads to vertex 
functions 



ax ■■■ajv;/3 



^({Ar}, {*},*;). (A6) 



For dimensions larger than 2 all perturbation contri- 
butions now have only singularities at least at d = 4. 
Further it may be convenient to introduce the correla- 
tion length instead of the temperature distance. In the 



present case two correlation lengths 



r H _ r M ! > = ^ In ({ Af }, {u}, fc 



,(A7) 



fc=0 



d 



q({Ar},{u}) = -^lnt^ \{Ar},{u},k 



(A8) 



k=0 



have to be introduced. Inserting the reversed equations 
Af_L = Ar ± ({£},{u}) and Af|| = Afy ({£}, {u}) into 
161) leads to functions 



r 



(A9) 



Introducing the correlation lengths in the vertex func- 
tions is a resummation procedure which removes the ex- 
panded contributions of the correlation length from the 
vertex functions. This leads to expressions which are 
simplified considerably. Moreover this is true for the cal- 
culations in dynamic models (see for instance part II). 
Within statics the two point functions reveal then the 
general structure 



f(2 
1 _L_L 



f (2,0) 



k) = 



(A10) 



(All) 



In two loop order the fc-independent function f± is 

2 



/x({f},{fi}) = l- 
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1|l,2„ T,±±\ 



-^<VvDZ^\{Z},0) (A12) 



where we have introduced the short notation 

dAfok) 



VfcS^O) 



dk 2 



(A13) 



fe=0 



The /c-dependent function g± reads in two loop order 

g ± ({Z},{u},k) = l 



nj_ + 2 



-^^(p^Hifrk) -Dt ±ll (m,0)) . (A14) 

In order to obtain /y and g\\ one only has to interchange 
_L and ||. The two loop integral D$ is defined as 



J 



D a r 3ak m,k) 



k> Jk» + fc' 2 )(^ 2 + fc" 2 )(& 2 + (k + k' + k") 2 ) 



(A15) 
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In the limit k — > we have 

A(t,k)-A(Z,Q) dA(£,k) 



lim 



k 2 



elk 2 

OK k=0 

V fe2 A(C,0) (A16) 



Applying the limit to g ai (see (|A12|) and (|A14|) ) one ob- 
tains 

Km ({£},{*}, A) = Lm.M) (A17) 

A: — >u 

and the two point vertex functions reduce in the asymp- 
totic region to 

(A18) 

Because the poles do not depend on k, the two functions 
f ai and g ai contain the same pole terms. Thus the two 
point functions each with may be renormalized by scalar 
renormalization factors which remove the poles from the 
functions f ai . 



The expansion coefficients a, in (|B2[) explicitly depend 
on couplings and coefficient Cyfc (|B1[) . Obviously, for 
t = 1 the function (|B1[) reproduces the initial RG func- 
tion ([B2]) : F(l) = /({it}). Then, the function dH2J> is 
resummed as a single variable function and further eval- 
uated at t = 1 to recover (|B1[) . To perform the resum- 
mation we use the Pade-Borel technique^ Namely, as- 
suming factorial growth of the expansion coefficients 
we define the Borel transform 3 - 3 - of (|B2[) by: 



F B (t) 



L 

£ 

i=0 



Oi/T(i + l)t l 



(B3) 



where T(x) is Euler gamma-function. Analytical con- 
tinuation of function (|B3[) is achieved by representing it 
in a form of a Pade approximant. 34 In our case, work- 
ing within a two-loop approximation we use the diagonal 
[1/1] Pade approximant: 



F B (t) ~ [!/!](*). 



(B4) 



APPENDIX B: RESUMMATION 

In this appendix we describe a procedure we use to 
resum divergent expansions for the two-loop RG func- 
tions. Starting from the RG function that has a form of 
truncated polynomial in renormalized couplings: 

L 

f({u})= J2 dikuiului (Bl) 
i,j,k=a 

one first represents it in a form of a resolvent series 3 - 1 * in 
variable t: 

L L 

F(t)= £ Cijfc«H«x* <+i+fc = Z)ai(M.{c})* < - 

(B2) 



Finally, the resummed function is obtained via an inverse 
Borel transform: 



F res = / [l/l](t)e _ *. (B5) 
Jo 



The procedure described above was used to analyze the 
RG flows and exponents. Note however, that the inverse 
Borel transform (|B5|) is well defined, when no poles in the 
denominator of Pade approximant (IB4|) appear. Other- 
wise one may estimate its principal value. The poles do 
not appear for a sign-alternating series (as the series for 
the /?- functions arc). To deal with sign alternating se- 
ries during an evaluation of critical exponents, we have 
resummed the functions 2 — £± . 
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